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Abstract A 

singularly perturbed linear system of second order partial differential equa- 
tions of parabolic reaction-diffusion type with given initial and boundary 
conditions is considered. The leading term of each equation is multiplied by 
a small positive parameter. These singular perturbation parameters are as- 
sumed to be distinct. The components of the solution exhibit overlapping 
layers. Shishkin piecewise-uniform meshes are introduced, which are used in 
conjunction with a classical finite difference discretisation, to construct a 
numerical method for solving this problem. It is proved that the numerical 
approximations obtained with this method are first order convergent in time 
and essentially second order convergent in the space variable uniformly with 
respect to all of the parameters. 



1 Introduction 



The following parabolic initial-boundary value problem is considered for a 
singularly perturbed linear system of second order differential equations 

du <9 2 u 

— E-—r + Au = f , on f2, u given on T, (f ) 

ot ax z 

where Q = {(a;, t) : < x < 1, < t < T}, J2 = Q U r, r = r L ur B ur R 
with u(0,i) = <p L (t) on r L = {(0,t) :0<t<T}, u(x,0) = 4> B {x) on r B = 
{(x,0) : < x < 1}, u(l,t) = cf> R (t) onr fl = {(l,t) : < t < T}. Here, for 
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all (x, t) G Q, u(x, t) and i(x, t) are column n — vectors, E and A(x,t) are 
n x n matrices, E = diag(e), £ = (ei, • • • , e n ) with < £i < 1 for all 
i = 1, . . . , n. The are assumed to be distinct and, for convenience, to have 
the ordering 

£!<••■<£„. 

Cases with some of the parameters coincident are not considered here. 
The problem can also be written in the operator form 

Lu = f on 17, u given on r, 

where the operator L is defined by 

L = d E° 2 I A 
dt dx 2 

For all (x, t) € O it is assumed that the components aij(x, t) of A(x, t) satisfy 
the inequalities 



aa(x,t) > / jajj(x,t)\ for 1 < i < n, and ay(x,t)<0 for i ^ j (2) 

3=1 

and, for some a, 

n 

< a < min (> dijix, t)). (3) 

l<i<n j = l 

It is also assumed, without loss of generality, that 

, \f0L 

m&x^-<v—. (4) 

l<i<n 

The reduced problem corresponding to ([lj is defined by 

+ Au = f, on fl, u = u on F B - (5) 

at 

The norms || V ||= maxi<fe<„ \Vk\ for any n- vector V, || y \\d— sup{|y(a;, t)\ : 
(x,t) £ D} for any scalar- valued function y and domain D, and || y || = 
maxi<fc<„ || yk \\ for any vector- valued function y are introduced. When 
D = Q or Q the subscript D is usually dropped. Throughout the paper C 
denotes a generic positive constant, which is independent of x, t and of all 
singular perturbation and discretization parameters. Furthermore, inequali- 
ties between vectors are understood in the componentwise sense. Whenever 
necessary the required smoothness of the problem data is assumed. 
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For a general introduction to parameter-uniform numerical methods for 
singular perturbation problems, see Q], [5] and [3]. The piecewise- uniform 
Shishkin meshes f] M ' N in the present paper have the elegant property that 
they reduce to uniform meshes when the parameters are not small. The prob- 
lem posed in the present paper is also considered in [5] , where parameter uni- 
form convergence is proved, which is first order in time and essentially first 
order in space. The meshes used there do not have the above typical property 
of Shishkin meshes. The main result of the present paper is well known in the 
scalar case, when n = 1. It is established in [3] for the case n = 2. The proof 
below of first order convergence in the time variable and essentially second 
order convergence in the space variable, for general n, draws heavily on the 
analogous result in [B] for a reaction-diffusion system. 

The plan of the paper is as follows. In the next two sections both stan- 
dard and novel bounds on the smooth and singular components of the exact 
solution are obtained. The sharp estimates for the singular component in 
Lemma [7] are proved by mathematical induction, while interesting orderings 
of the points established in Lemma El In Section 4 piecewise- uniform 

Shishkin meshes are introduced, in Section 5 the discrete problem is defined 
and the discrete maximum principle, the discrete stability properties and a 
comparison result are established. In Section 6 an expression for the local 
truncation error is derived and standard estimates are stated. In Section 7 
parameter-uniform estimates for the local truncation error of the smooth 
and singular components are obtained in a sequence of theorems. The sec- 
tion culminates with the statement and proof of the essentially second order 
parameter-uniform error estimate. 



2 Standard analytical results 

The operator L satisfies the following maximum principle 

Lemma 1. Let A(x,t) satisfy (dP and Let ib be any function in the 
domain of L such that ij) > on T. Then Lxp(x) > on fl implies that 
ip(x) > on Q. 

Proof. Let i*, x* , t* be such that ipp (x* , t*) = min^ mhtj)ipi(x, t) and assume 
that the lemma is false. Then ipi» (x* , t*) < . From the hypotheses we have 

(x*,t*) £ r and ^^(x*,t*) > 0. Thus 

(LiP(x*,t*)) t , = ^l( x *,f)- ei .^(x* J t*)+^o i . j3 -(a : * 1 i*)^(a : *,i*) < 0, 
which contradicts the assumption and proves the result for L. ■ 
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Let A(x, t) be any principal sub-matrix of A{x, t) and L the corresponding 
operator. To see that any L satisfies the same maximum principle as L, it 
suffices to observe that the elements of A(x,t) satisfy a fortiori the same 
inequalities as those of A(x,t). 

Lemma 2. Let A(x,t) satisfy H)) and |3J). Ifxp is any function in the domain 
of L, then, for each i, I < i < n and {x,t) £ £2, 

\ipi(x,t)\ < max 1 1| V ll-r, ^ II Li/> ||| . 

Proof. Define the two functions 

9 ± (x,t) = max||| t/> \\r, — || Lt/> |||e ± ip(x,t) 

where e = (f , . . . , 1) T is the unit column vector. Using the properties 
of A it is not hard to verify that 6 ± > on r and L0 ± > on ft. It 
follows from Lemma [T] that 6 ± > on I? as required. ■ 

A standard estimate of the exact solution and its derivatives is contained in 
the following lemma. 

Lemma 3. Let A(x,t) satisfy (0) and f3J) and let u be the exact solution of 
(QP- Then, for all (x, t) £ Q and each i = 1, . . . , n, 



\^(x,t)\ <C(||u||r + EUll^H)' * = 0,1,2 



\^r(x,t)\ <CeT(||u|| r + ||f|| + ||Il|), 1 = 1,2 



1-2 



\%£?(x,t)\ < Ce-'e— (\\u\\ r + ||f|| + \\§ \\ + \\^\\ + \\§^\\), I 

\^m(x,t)\ < Ce^(\\u\\ r + ||f|| + Hfll + II0H), I = 2,3. 

Proof. The bound on u is an immediate consequence of Lemma [2j Differenti- 
ating (JIJ partially with respect to once and twice, and applying Lemma [2] 
the bounds and 0^ are obtained. To bound ^f, for all i and any (x,t), 
consider an interval I = (a, a + y/Sj) such that x £ /. 
Then for some y £ L and t £ (0, T) 

dm, . Ui(a + Js~i,t) - Ui(a,t) 

&>' ) = ^ 



|^(y,t)|<C £ ^||u||. (6) 



Then for any x £ L 
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rx pa 



dm , , dui , . f x d Ui(s,t) , 

^(*.*) = ^(y,*) + y B 

i^ ( , )t )i<i^ ( ^)i + ^r ( || U || r+ i| f || + ii|ii) rf , 



Using dHJ) in the above equation 

|^(x,OI<C e ^(||u|| r + ||f|| + |||| 
Rearranging the terms in ([T]), it is easy to get 

' i>2 " i '<Ce7 1 (\\u\\r + \m\ + \&\). 



1 dx 2 1 - 1 v " IM at 1 

Analogous steps are used to get the rest of the estimates. ■ 

The Shishkin decomposition of the exact solution u of([T])is u = v + w 
where the smooth component v is the solution of 

Lv = f in J?, v = Uo on r (7) 

and the singular component w is the solution of 

Lw = in ]7, w = u — v on r. (8) 

For convenience the left and right boundary layers of w are separated using 
the further decomposition w = w L + where Lw L = on 17, w L = w 
on Tl, w l = on Tb U -Tr and Lw R = on J?, xv R = w on w R = 

on r L ur B . 

Bounds on the smooth component and its derivatives are contained in 

Lemma 4. Let A(x,t) satisfy 0j and Then the smooth component v 
and its derivatives satisfy, for all (x, t) € Q and each i — 1, . . . , n, 

\^r{x,t)\ < C for 1 = 0,1,2 

\f^(x,t)\ < C(l + £ H) for Z = 0,1,2,3, 4 

\ 1 *% m (x,t)\<C for 1 = 2,3. 

Proof. The bound on v is an immediate consequence of the defining equations 
for v and Lemma ((2]). Differentiating the equation ([7]) twice partially with 
respect to x and applying Lemma [2] for , we get 
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i£i*°»+ii£ii>- M 

Let 

— — (x*,t*) = ||-=— II for some i — i* , x = x* , t = t* . (10) 
ox ox 

Using Taylor expansion, it follows that, for some y € [0, 1 — x*\ and some 
?7 S (ir*,2;* +y) 

+ y,t*) = v l ^x*,t*)+y^(x*,t*) + ^^( v ,t*). (11) 



Rearranging (fTTj) yields 



g«g / * ,*s _ ^»Qc* +V,t*) - Vj* (x*,t*) _ y d 2 Vj, 
dx [X ' >- y 2 fe 2 W ' ' 

l ^ (l .,, )l£ a l|v|l + fll g ll . (12) 



Using UTDJl and (pj) in ([9]l. 



This leads to 



or 



l^l<C(l + -||v|| + -||^| 



(i-?)ii£ii^ a+ > 



»ft?^ c - (13) 

Using JT3]) in (HU) yields 

ll?ll<C 
ox 

Repeating the above steps with it is easy to get the required bounds on 

<9 3 v d v 

the mixed derivatives. The bounds on -^—5-, -7— 7 are derived by a similar 

ax J ax 4 

argument. ■ 



3 Improved estimates 



The layer functions Bf , Bf 1 , Bi, i = 1, . . . , n, , associated with the solution 
u, are defined on [0, 1] by 
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B L (x) = e - X ^J7 % B R {x) = B L {1 _ xl Bi[x) = B L (x) + B R {x y 

The following elementary properties of these layer functions, for all 1 < i < 
j < n and < x < y < 1, should be noted: 

(a) flftaO < B^x), Bf{x) > B?(y), < Bf{x) < 1. 

(b) Bf- (x) < Bf(x), Bf-(x) < B?(y), < B?(x) < 1. 

(c) Bi(x) is monotone decreasing (increasing) for increasing x E [0, ^]([^, 1]). 

(d) Bi(x) < 2Bf{x) for x G [0, §], B^x) < 2Sf (as) for a; e [±, 1]. 

Definition 1. For £?f, each i,j, 1 < i ^ j < n and each s, s > 0, the 

f s) 

point ^ is defined by 



(14) 



It is remarked that 



£f(l-zg) ffil-sj'?) 



(15) 



» 

HJUCUCOO ui tile puiiito ./* 

Various properties are also established. 



In the next lemma the existence and uniqueness of the points xfl are shown. 



Lemma 5. For all such that 1 < i < j < n and < s < 3/2, the points 
are uniquely defined and satisfy the following inequalities 

Moreover 

x ( /) < x$ 1;j , if i+Kj and &JJ < if i < j. (17) 



< 2s ^ and a>> € (0, \) if i < .,. 



(18) 



Analogous results hold for the Bf , B^ and the points 1 — . 



Proof. Existence, uniqueness and (|16[) follow from the observation that the 
ratio of the two sides of (jl~4l , namely 



Bt(x) ej e s 3 1 l 

— L - ! ^- — r^r = — exp (— Vaa;(— ^= )) 

el Bf{x) el PK V ^ ^e]" 
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is monotonically decreasing from the value p- > 1 as x increases from 0. 

(s) . . . 

The point x\ is the unique point x at which this ratio has the value '. 
Rearranging (|14[) , and using the inequality In x < x — 1 for all x > 1 , gives 

2 S ln(^E) [7 



7=) 



' <2*«p, (19) 



which is the first part of (|T9|) . The second part follows immediately from this 
and ®. 

To prove (JT7J) , writing = exp(— pk), for some pfc > and all k, it follows 
that 

(s) _ 2.s(p. t -pj) 



JJ Va(exppi - exppj) 



The inequality x\ i < • is equivalent to 



Pi - Pj Pi+i ~ Pj 



exppi - exppj exppi + i - exp p., 

which can be written in the form 

(Pi+i ~ Pj) exp(pi -pj) + {pi - Pi+x) - (Pi - Pj) exp(p 4+ i - pj) > 0. 

With a = Pi — Pj and b = Pi+\ — Pj it is not hard to see that a > b > and 
a — b = pi — Pi+i- Moreover, the previous inequality is then equivalent to 

exp a — 1 exp 6—1 
a b 

which is true because a > b and proves the first part of (flTj) . The second part 
is proved by a similar argument. 

The analogous results for the B^, B^- and the points 1 — x\j are proved by 
a similar argument. ■ 

In the following lemma sharper estimates of the smooth component are pre- 
sented. 

Lemma 6. Let A(x, t) satisfy (0) and f3J). Then the smooth component v 
of the solution u of ^ satisfies for all i = 1, • • • , n and all (x,t) e Q 

l§(M)|<c(l + £^| 6.1 = 0,1,2,3. 
\ Q=i £ q 

Proof. Define two barrier functions 

rp ± (x,t) = C[l + B n (x)]e± -^-(x,t), i = 0,1,2 and (x, t) G J2. 
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We find that, for a proper choice of C, 

i,f(0,t) = C ± ^(0,t) = C>0 
V>f(M) = C ± ^l(i,t) = c>0 

t/>t(x,0) = C[l + B n (x)} ± d ^ X) =C[l + B n (x)} ± C> 

as <p b (x) £ C( 2 )(r 6 ) and (Ltp^^t) > 0. 
Using Lemma [1] we conclude that 

l-^jMI <C[l+B n {x)\ for Z = 0,1,2. (20) 
Consider the equation 

1 [ d X 2 )h ~ dx* dx dx dx* 3 [ > 

with 

0,0,0. 0,0,1,0 = o.0(,o, = m 

For convenience, let p denote ^pr- Then 

Lp = g with p(0,i) = 0, p(l,i) = 0, p(x,0)=s (23) 

where 

- _ o ^J-i _ ^il v . an d s . - (g) 

1 9a; 2 9a; 9a; ^— ^ 9a; 2 J 9a; 2 

Let q and r be the smooth and singular components of p given by 

Lq = g with q(0,t) = p o (0,t), q(l,t)=po(l,t) ) q(x, 0) = p(x, 0) (24) 
where po is the solution of the reduced problem 

+ Ap = g with p (x, 0) = p(x, 0) = s. 

Now, 

Lr = 0, with r(0, t) = -q(0, t), r(l, i) = -q(l, t), r(x, 0) = 0. (25) 
Using Lemma |4] and Lemma O we have for i = l,...,n and (x,i) G i? 
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lg(M)|<C 

and 



dx y/Fi 
Hence, for (x, t) € Q and i = 1, . . . , n, 

&\ = \ d f\<C[l + B p + - + B p]. (26) 
dx A dx y/Ei y/e^ 

From and (gU), we find that for I = 0, 1, 2, 3 and (x, t) G 77 

l^rl < C[l + e^B.ix) + ■■■+ e]-*B n {x)]. 

Remark : ft is interesting to note that the above estimate reduces to the 
estimate of the smooth component of the solution of the scalar problem given 
in [I] when n = I. 

Bounds on the singular components w L , of u and their derivatives are 
contained in 



Lemma 7. Let A(x,t) satisfy Hp and f3J). Then there exists a constant C, 
such that, for each (x, t) £ fi and i = I , . . . , n, 

\^f(x,t)\<CB^(x), for Z = 0,l,2. 
I#(M)I<C£^^R for 1 = 1,2. 

Analogous results hold for wf" and its derivatives. 

Proof. To obtain the bound of w L , define the functions ^^(Xjt) = 
Ce at B^(x) ± wf{x,t), for each i = l,...,n. It is clear that ^^(0,*), 
ipi^fa, 0), ^^(l, t) and (Lxp ± )i(x, t) are non-negative. By Lemma 1, ip^ix, t) > 
0. It follows that \wf\ < Ce at B^{x) 
or 

\wf\< CB%(x). (27) 
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To obtain the bound for 5^- , define the two functions 8i ± (x,t) — CB^ {x) ± 

-^-{x,t) for each i = l,...,n. Differentiating the homogeneous equation 
satisfied by wf 1 , partially with respect to t, and rearranging yields 



at 2 ^ l dx 2 dt + ^ lJ at ~ at j ' 

3 = 1 

and we get 

\L-^\<CB n L (x) 

l^(o,*)l<l^(o,*)l + l^(o )t )| = |^|<c 

l^(l,*)l<l^(l,*)l + l^(M)l = l^l<C 

l^,o)| <|^M| = o. 

By Lemma [5J it follows that 

\^\<CB n L (x). (28) 
Now the bound for dx £ t is obtained by using Lemma ([3]) and Lemma (U) 

, a 2 w 1 t , , a 2 u,, , , a 2 vi , 



<9x<9i dxdt 
U 2 >c<: „ -i „,, ,,<9f,, ,,5 2 f, 



<C£jT u r + f + - +1 



Similarly, 



dxdt 1 " 1 v " l|J dt dt 2 



1^1 ^aMir + i m + + llgl 



The bounds on , Z = 1,2,3,4 and i = 1, . . . , n are derived by the method 

of induction on n. It is assumed that the bounds ^-J-, , and a 
hold for all systems up to n — 1. Define w L = then w L 

satisfies the system 

dw L ~ d 2 w L ~ . r 



at dx 2 

with 



w L (0,t) = u(0,t) - u o (0,i),w L (l,t) 



12 
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w l (x,0) = u(x,0) - uo(x,0) = 4>b(x) - 4>b{x) 



0. 



Here, E and A are the matrices obtained by deleting the last row and column 
from E, A respectively, the components of g are gi — —ai tn w^ for 1 < i < 
n — 1 and iio is the solution of the reduced problem. Now decompose w L 
into smooth and singular components to get w L = q + r, = + ff ■ 

By induction, the bounds on the derivatives of ~w L hold. That is for i = 
l,...,n- 1 



dw L 
I _l I 

1 dx 1 

d 2 w L 
I i I 

1 dx 2 1 
d 3 w L 

I _ 1X1 % I 

1 dx 3 1 

<9V L 



a^ 4 



(29) 



Rearranging the n equation of the system satisfied by u>„ yields 



dx 2 



Of 



Using (HID and (28J) gives 



1 &B 2 



< Ce-'B^x). 



(30) 



Applying the mean value theorem to at some y, a < y < a + 



9W "-(y t) = W "^ + ~ 



<9.c 



Using (f27|) gives 



r)m L C 

-.(y,t)\<—(B^a + ^ n ) + B^(a)). 



dx 



So 



Again 



, dw, 



C 



(31) 



"flf ( a: '*) = -ftf(»' < ) + (f- a; )-^f('/.*). 2/<^<*- (32) 
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Using flU]) and ([ST]) in ([32]) yields 

\?g(x,t)\<C[e?BZ(a)+ E ?BZ(T,)] 

— r--~ R L (t\ B n (°) 

— Ot„ JD n \d,) gl ^ 

= Ce7~B%(x)e^nVZ/Ve n , 



Therefore 



,dwi 



(x,t)\ < Ce n 2 B^(x). (33) 
ox 

Now, differentiating the equation satisfied by partially with respect to x, 
and rearranging, gives 

£ " 9x3 " ^ + 2^ a »« ^ + a™ ^ +2^ 5a; % ■ 

q=l q=l 

The bounds on and then give 

Q3.„L 



-^\<C^e q ^B^x). 

9=1 



Similarly 

9=1 



' 9a; ' "SS 5 " ' ~d^~ aild "35^ 

C £ n 1 2g=i ~~7^" respectively. Introducing the functions ^ ± (x, t) — CB^ (x)e ± 

q(ir,i), it is easy to see that i/? ± (0,i) = Ce ± q(0, t) > 0, i/> (l,t) = 
CB%(l)e ± > 0, '0 ± (^,O) = CB%(x)e ± > and 



Using the bounds on u>„ , , , and it is seen that 

3g d*g 3fg KniinHorl Hir HR^M V™ ^_fl 



W ± ) j (x,t) = fJ(-£ i -+ra ii )B n i (^) ± CB^(x)>0. 



Applying Lemma 1, it follows that ||q(a;,t)|| < CB^(x). Defining barrier 

functions 8 ± (x,t) = Cen B^(x)e ± ^ and using Lemma 3 for the problem 

satisfied by q, the bound required for ^ and is obtained. By induction, 
the following bounds for r are obtained for i = 1 , . . . , n — 1, 
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\ d2ri \ <c\ B i {x) i i g "- l(a:) ], 

d 3 r ■ Bt(x) B^jx) 

\-Q^J\ ^ M 1 1 1 I J' 



1 c n-l 



Combining the bounds for the derivatives of qt and it follows that, for 
i = 1, 2, . . . , n 



I dx' I — I a^r I + I f^r I 

1^1 <CE n H ^ for 1 = 1,2 



and | e ^|<c E: ££(fl. 

Recalling the bounds on the derivatives of completes the proof of the 
lemma for the system of n equations. 

A similar proof of the analogous results for the right boundary layer functions 
holds. ■ 



4 The Shishkin mesh 

A piecewise uniform Shishkin mesh with M x N mesh-intervals is now con- 
structed. Let CI? = {t k }f =l , Q% = {xj}^ 1 , Jlf = {tfc}^, < = 

{x,}f =0 , n M > N = x n M ' N = 77f x q n x and r M ~ N = rnQ M ' N . 

M 

The mesh fi t is chosen to be a uniform mesh with M mesh-intervals on 

[0, T]. The mesh Q x is a piecewise-uniform mesh on [0, 1] obtained by divid- 
ing [0, 1] into 2n + 1 mesh-intervals as follows 

[0, <7i] U • • • U K_i,<T„] U (cr„, 1 - a n \ U (1 - <r„, 1 - <r„_i] U ■ ■ ■ U (1 - (71, 1]. 

The n parameters oy, which determine the points separating the uniform 
meshes, are defined by 
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oy, = min <{ -, 2 J — InN \ (34) 



4' V a 



mini^-,2J^\nN}. (35) 
' 2 V a 



and for r = 1 , ... , n — 1 



Clearly 

< a 1 < ... < o n < ^ < 1 - o n < ... < 1 - oi < 1. 

Then, on the sub- interval (o n , 1 — o n ] a uniform mesh with mesh- 
intervals is placed, on each of the sub-intervals (o ri o r+ i] and (1 — oy+i, 1 — 
ay], r = l,...,n — 1, a uniform mesh of 2 „- r +% mesh-intervals is placed 
and on both of the sub-intervals [0, o\] and (1 — oi, 1] a uniform mesh of 

N 
W+ 1 



A mesh-intervals is placed. In practice it is convenient to take 



N = 2 n+p+1 (36) 

for some natural number p. It follows that in the sub-interval [oy-i, oy] there 
are N/2 n ~ r+3 = 2 r+p ~ 2 mesh-intervals. This construction leads to a class of 
2 n piecewise uniform Shishkin meshes f2 M,N . Note that these meshes are not 
the same as those constructed in [5] 

The following notation is introduced: hj — Xj — Xj—i, J — {xj : D + hj = 
hj+i — hj 7^ 0}. Clearly, J is the set of points at which the mesh-size changes. 
Let R — {r : oy 6 J}. From the above construction it follows that J is a 
subset of the set of transition points {oy}" =1 U {1 — <r r }" =1 . It is not hard to 
see that for each point Xj in the mesh-interval (oy_i, oy], 

hj = 2 n - r+3 N- 1 {o r - o r -i) (37) 

and so the change in the mesh-size at the point oy is 

D+h r = 2 n - r+3 (d r - d r -x), (38) 

where d r = <T ' 2 +1 — oy for 1 < r < n, with the conventions do — 0, o n+ \ = 1/2. 
Notice that d r > 0, that Q M ' N is a classical uniform mesh when d r = for 
all r = 1 ... n and, from (|38[) . that 

D+h r < if d r = 0. (39) 



Furthermore 

o r < C \fs~r In N, l<r<n, (40) 

and, using ([37]), gQl), 
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Also 

ay = 2~ (s ~' r+1) cr s+ i when d r = ■ ■ ■ = d s = 0, 1 < r < s < n. (42) 

The geometrical results in the following lemma are used later. 

Lemma 8. Assume that d r > and let < s < 2. TTien the following 
inequalities hold 

B^(a r )^N- 2 . (43) 

4-l,r < °V - /V for 1 < r < n. (44) 

•Si-(cTr) 1 

g s < — for 1 < 9 < n, 1 < r < n. (45) 
^(oy - h r ) < CB^(a r ) for 1 < r < q < n. (46) 



Proof. The proof of (|43[) follows immediately from the dchnition of ay and 

the assumption that d r > 0. 

To verify (|44)) note that, by Lemma [5] and (|36|). 

(s) „ j£r _ so r _ say a r 

Xr ~ hr < S V 77 ~ InTv ~ (n+p+l)ln2 ~ T' 

Also, by (1551) and (|37|) . 

/ V = 2"- r+3 A^- 1 (a r - oy^) = 2 2 - r -f( ( r r - oy^) < ° r ~ 2 f7r ' 1 < ^. 
It follows that x r _ x r + /i r < ay as required. 



To verify (|45j) note that if g > r the result is trivial. On the other hand, if 
q < r, by (|31)l and Lemma [5l 

iff M < Bj-(arffi) _ g r L (x^) < 1 

Eq £g S r € r 

Finally, to verify (14T>1) note, from (j3"T)l . that 

ft r = 2™- r+3 A- 1 (cr r - ay_ x ) < 2"- r+3 iV-V r = 2 n - r+4 J—N- 1 InJV. 



But 



so 
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'— h r < J— 2 n - r+4 ^- 1 lniV < 2 n - r+4 N- 1 \nN < C, 
since r < q. It follows that 



as required. 



5 The discrete problem 

In this section a classical finite difference operator with an appropriate 
Shishkin mesh is used to construct a numerical method for ([T]), which is 
shown later to be essentially second order parameter-uniform. It is assumed 
henceforth that the problem data satisfy whatever smoothness conditions are 
required. 

The discrete initial-boundary value problem is now defined on any mesh by 
the finite difference method 

DtU - E%U + AU = f on Q M ' N , U = u on r M ' N . (47) 

This is used to compute numerical approximations to the exact solution of 
([1}. Note that (|4T|) . can also be written in the operator form 

L M ^U = f on Q M > N , U = u on r M > N , 

where 

L ALN = Dt -E81+A 
and , 8%, and D~ are the difference operators 



S 2 x U(x Jltk ) 



U(xj,t k ) - U(xj,t k -i) 
tk — tk-i 

Dj\J{x v t k )-D-\J{x ,t k ) 



(x J+1 -Xj-x)/2 



-f TT( ... _ U(x j+1 ,t k ) -V(xj,t k ) 



x o x j-t 

The following discrete results are analogous to those for the continuous case. 
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Lemma 9. Let A(x,t) satisfy |J|) and (0). Then, for any mesh function \&, 
the inequalities * > on T M ' N and L M ^* > on Q M > N imply that 

—M,N 

* > 071 Q 

Proof. Let i*,j*,k* be such that (xj* , t k * ) = min^ min^fe &i(xj t t k ) and 
assume that the lemma is false. Then <Zv (xj* ,t k *) < . From the hypotheses 
we have f ^ 0, N and fy. (xj, , t k , ) - fy, (a; J -»_i, i fe ») < 0, - 
(xj. , > 0, so 5lWi^(x^,t k ,) > 0. It follows that 

n 

(L^*^*,^.))^ = -e^o^w^j^tfcO+X! "' q(xj*,t k *)\P q (Xj*,t k *) < 0, 

9=1 

which is a contradiction, as required. ■ 

An immediate consequence of this is the following discrete stability result. 

Lemma 10. Let A(x,t) satisfy 0/ 1 and f3J). Then, for any mesh function \& 

071 ft, 

|| ¥{ Xj ,t k ) || < max(||*|| r M, N ,i||L M ^*| 
Proof. Define the two functions 

& ± (x j ,t k ) = max{||*|| r M,», i||L M ' N *||}e± 9(x jt t k ) 

a 

where e = (1, ... , 1) is the unit vector. Using the properties of A it is not 
hard to verify that ± > on r M > N and L M > JV ± > on Q M ' N . It follows 
from Lemma IH1 that ®^ > on Q ' . ■ 

The following comparison result will be used in the proof of the error estimate. 

Lemma 11. (Comparison Principle) Assume that, for each i = 1, ... ,n, 
the mesh functions <& and Z satisfy 

\Zi\<$i, on r M > N and |(L MJV Z),;| < (L M ^*) 4 on f2 M > N . 

Then, for each i = 1, ... , n, 

\Zi\<$i. 

Proof. Define the two mesh functions by 

* ± = $±Z. 



Then, for each i = 1, ... , n, satisfies 
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(3^)i>Q, on r M ' N and \(L M ' N Z) i \ < (L ALN &)i on f2 M ' N . 
The result follows from an application of Lemma [9] ■ 

6 The local truncation error 

From Lemma [TOl it is seen that in order to bound the error ||U — u|| it suffices 
to bound h M ' N (U — u). But this expression satisfies 

L M - N (U - u) = L M ^(U) - L M < w (u) = 

f - L M > N (u) = L(u) - L M ^(u) = (L - L M > N )u. 

It follows that 

L^(U-u) = (|- J D t -)u-£(^-^)u. 

Let V,W L , ~W R be the discrete analogues of v,w L ,w R respectively. Then, 
similarly, 

L^(V-v) = (|-i5 t -)v-^-^)v ) 

l/^(W* - w*) = (| - A")w fl ' - - £)w« 

and so, for each i = 1, ... , n, 

|(L A/ ^(V - v)),| < |(| - D t ->| + - (48) 

|(L M ^(W L - w L )) 4 | < |(A _ D t ~)«;f | + \ £l (^ - 5l)wfl (49) 

|(L A ^(W* - w»)) 4 | < |(| - D t -)tuf | + | £ ,(^ - (50) 

Thus, the smooth and singular components of the local truncation error can 
be treated separately. Note that, for any smooth function ip, the following 
distinct estimates of the local truncation error hold: 
for eacb(Xj,t k ) £ O l ' m 

\(i-Dt)if(xj,t k )\ < C{t k -t k - X ) max \^t( Xj ,s)\, (51) 
oi s e [tfc_i,tfc] ox 
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and 

1(^2 -Wfe-.**)! < C(a; J+ i -^_!) max 1^-3(5, t fc )|. (53) 
Assuming, furthermore, that Xj £ J, then 

d 2 d A 4> 
K^2 -W(a;j,*fc)l < ^(^'+1 -^-i) 2 max| — (s,t fc )|. (54) 

Here Jj = [xj-i, 



7 Error estimate 

The proof of the error estimate is broken into two parts. In the first a theorem 
concerning the smooth part of the error is proved. Then the singular part of 
the error is considered. A barrier function is now constructed, which is used 
in both parts of the proof. 

For each r € R, introduce a piecewise linear polynomial 9 r on f2, defined by 

X 

— , < x < a r . 

(T r 

e r (x) = { 1, O r < X < 1 - <7 r . 
' 1 — X 

, 1 — <7 r < X < 1. 

O r 

It is not hard to verify that 

{a6 r (xj), if Xj ^ J 
<j r (h r + flr+l) 

On the Shishkin mesh S7 M < N define the barrier function 3? by 

&(xj,t k ) = CiM- 1 + N- 2 + N- 2 (\nN) 3 J2 6 r(x 3 )]e, (56) 

r£-R 

where C is any sufficiently large constant. 
Then, on Q M > N , * satisfies 

< <&i(xj, tk) < C(M _1 + Af _2 (ln N) 3 ), 1 < i < n. (57) 

Also, for Xj ^ J, 
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(L M > N <f>)i( Xj ,t k ) > C'iM- 1 +N- 2 (lnN) 3 ) (58) 
and, for Xj £ J, using (|4T|) . (|55)l . 

{ C{M~ l +N- 2 + - F M=N- 1 biN), if D+h r > 0, 

(59) 

The following theorem gives the required error estimate for the smooth com- 
ponent. 

Theorem 1. Let A(x, t) satisfy Hp and (0|). Let v denote the smooth com- 
ponent of the exact solution from (OP and V the smooth component of the 
discrete solution from |77p . Then 

||V- v|| < C'iM- 1 + N- 2 (hiN) 3 ). (60) 

Proof. It suffices to show that 



\(L M > N $)i(xj,t k )\ 



< C, (61) 



for each i = 1, ... , n, because an application of the Comparison Principle 

then yields the required result. 

For each mesh point Xj either Xj ^ J or Xj £ J. 

Suppose first that Xj ^ J. Then, from (j58")l , 

(L M ' N ^( Xj ,t k ))i > CiM-'+N- 2 ) (62) 

and from (fSTj) . (|54|) and Lemma [4] 

((L^V - v)) t (x 3 ,t k )\ < C(t k - t k -i + (x j+ i - x^) 2 ) 

< CiM- 1 + (hj + h 3+l ) 2 ) (63) 

< CiM- 1 +N- 2 ). 

Then (|tTTj) follows from (|()2"j) and (j6"3")l as required. 

On the other hand, when Xj £ J, by (ISTl) . (|53|) and Lemma [6] 

|(i M ^(V - v))^-,**)! < C[M~ X + edhr + K +1 )(l + B ^~ K) )]. 



q=l 

(64) 

The cases i > r and i < r are treated separately. 
Suppose first that i > r, then it is not hard to see that 

\(L M > N {V-vM Xj ,t k )\ 

< C[M~ X + e l {h r + h r+1 )(l + -±=)] (65) 

< C[M~ l + (h r + h r+1 )^/el]. 
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Combining (|5^|) and , (|61l) follows using (jUJ) and the ordering of the Si . 
On the other hand, if i < r then £j < e r _i < e r . Also, either d r > or 
d r = 0. 

First, suppose that d r > 0. Then, by Lemma [51 



h r > Xqf r ; for i < q < r — 1 



and so, by Lemma [5] 



r-i 



y-^ B q (a r — h r ) ^ £,B r {r>, h, ) 



Combining this with (jMJ) gives 



v 7 ^ 



| (i M,JV (v _ v))^^.,^)! < dM' 1 + ^(h r + hr+i)]. (66) 



Combining (15^)) and (j6"6")l . (|61l) follows using (j^TJ) and the ordering of the S{. 
Secondly, suppose that d r — 0. Then d r _i > and D + h r < 0. Then, by 
Lemma [8j with r — 1 instead of r, 

<T r — h r > a r -\ > oy_i — /i r _x > x q r~i f° r * — 9 — r ~ ^ 
and so, by Lemma [5] 

^ £q(oy - hr) < c , -Br_i(qV ~ K) < g gr_i(Or-l) _ g N' 2 
q=i V^Q y/^r-l — 1 y/^r^l 



Combining this with (|64| and (|4TJ) gives 
|(LM,iv (v _ v)) . (a; . )tfc) | 

<C[M- 1 +£ 4 (/ lr + / lr+1 )(-g= + -^)] (67) 
< CfA-f- 1 + y /£~e;N- 3 In 7V + £ 4 7V-i In AT]. 

Combining (l59t and (|67)) , (j6"Tj) follows using the ordering of the e\ and noting 
that in this case the middle term in the denominator is used to bound the 
middle term in the numerator. 



Before the singular part of the error is estimated the following lemmas are 
established. 



Lemma 12. Assume that Xj ^ J. Let A(x,t) satisfy £||) and f3J). Then, 
Q M - N , for each 1 < i < n, the following estimates hold 



on 



\(L M ' N (W L - w L )M^,i fe )| < C(M~ 1 + {Xj+1 Xj - l)2 ). (68) 
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An analogous result holds for the wf~. 

Proof. Since Xj ^ J, from (|54|) and Lemma [3 it follows that 

| (jL m, N(w l _ w i Mxj)tk )\ = |((JL _ D -) _ _ «52))wf (^.tfc)! 

™ B L (s) 

< C(M : + - x^-i) 2 max ^ ' ' < 



< C(M _1 + 

as required. ■ 

The following decompositions are introduced 



r+1 



E 



W, L = > Wi m . 



m—1 



where the components are defined by 



I wr otherwise 



and for each m, r > m > 2, 

"pP on [0,x« 1>m ) 

r+1 

— otherwise 

r?=m+l 

and 

r+1 

w lA = wf - on [0, 1]. 

g=2 

Here the polynomials , for s = 3/2 and s = 1, are defined by 
and 



,1). , (1) N (x 

pi ^ = E -^irKrW r , • 

9=0 



Lemma 13. Assume that d r > 0. Let A(x,t) satisfy and J3|). Then, for 
each 1 < i < n, there exists a decomposition 
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r+1 



q=\ 



for which the following estimates hold for each q and r, 1 < q < r, 
1^(0^)1 <Cmin{l j}B^{ Xj ), 

q i 

< (> '" k ' q -( Xj ,t k )\ <Cmin{4^-^=K(2^ 



I ^ 3 (xj,t fc )l <Cmm{ ^ —375-, 2^ 

q=r+l £ q q=r+l 'V 9 



\^-('jM<c ^ — — 



q=r+l 



% Cr q 



Analogous results hold for the wf' and their derivatives. 

Proof. First consider the decomposition corresponding s = 3/2. 

From the above definitions it follows that, for each to, 1 < to < r, 

».,m =0 on [x^ +1 ,l]. 

To establish the bounds on the third derivatives it is seen that: for x € 
[x^Z+i ' 1] 1 Lemma [7j and x > /+\ i m Ply that 

q=l £ q q=r+l £ q 

for x £ [0, a^r+i]i Lemma [7] and x < xfl+\ imply that 



d 3 Wj, r+1 _ dVf (3/2) 



_3/2 — ^.3/2 — __3/2 

-9 



and for each to = r, ... , 2, it follows that 

for x E [x%™ +1 , 1], wl 3 l = 0; 

for x € [a^m-i m' i2m|ii Lemma [7] implies that 



Numerical solution of a parabolic reaction-diffusion system 25 

l-feS-(».*)l<l-^-(«.*)l+ E l^r^l 

q— rn+1 



3/2 - ~ 3/2 

q=l t 9 fc »i 

„(3/2) I 



for x £ [0,x) m _ 1 m ], Lemma [7] and x < x m _{ m imply that 

I g x 3 ( x »*)l = I Q x 3 ( x m-l,m)*)l 

n B L (r (3/2) \ B L (t {3/2) ) RLM 

^ q~ST^ q \ m—\,rrh) ^ Q JJ m\ Jj m—l,mi ^ Q D m\ x ) , 
— 3/2 — 3/2 — 



3/2 — ^3/2 — ^.3/2 

(3/2) nl a 3 ™,,! 



g__l £q £m £Y, 



farxe^M], ^ = 0; 
for x € [0, x\ 2 ], Lemma [3 implies that 

< 1^)1 + £ ,^«M), £ CE ^ < c^i. 

<j=2 g=l e <? £ 1 

For the bounds on the second derivatives note that, for each m, 1 < m < r 

r -- r (3/2) i] d 2 Wi,m o 

tor a; t Fm,m+i> L U Ox 3 — u > 

for x e [0,4 3 S +1 ], %-( S ,t)d S = 

C? 2 «'i m f (3/2) ,\ 9^t£i_m/ .\ d 2 Wj m f ,\ 

dx 2 \ X m-l,m' l ) — dx 2 \ X -> 1 ) ~ ~ dx 2 \ X i l ) 

and so 

(3/2) (3/2) 

This completes the proof of the estimates for s = 3/2. 

For the estimates in the case s — 1 consider the decomposition 

r+l 



E 



From the above definitions it follows that, for each m, 1 < m < r, 
on [x^n+ul]. 

To establish the bounds on the fourth derivatives it is seen that: 
for x £ [^ r . r -|-i, 1], Lemma [7] and x > x r r+1 imply that 
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q—1 Q q— r+1 ^ 

for x € [0,x^ +1 ], Lemma [7] and x < x^ +1 imply that 



™ B L (r t ' 1 ^ ) ™ B L (r W \ n B l (t) 

q—1 " q— r+1 " g=r+l ^ 

and for each m = r, ... , 2, it follows that 

„(i) il 9^,, 



for X € [ x m ,m+l' 1]' da; 4 — ^ 



for x € m , ^m^m+ili Lemma [7] implies that 



9=1 * 

for x € [0, a;^_i m ], Lemma [7] and x < m imply that 

\ £i g x 4 ( x >*)l = l £i g 4 ( x m-l,m)*)l 
" nLf (!) \ Rid 1 ) ^ uLr^ 

for x € [x^,l], = 0; 

for x € [0, x^], Lemma [7] implies that 

| ei ^ M )|<h^(M,l + Ek,^(M)| 

< c y5£M< c £fM. 

For the bounds on the second and third derivatives note that, for each m, 
1 < to < r : 

r r (1) -i 1 d 2 Wi_ m p. d 3 Wi ;m 

ior x t p mjm+ i, rj, — Q^t— — u — a a 3 ; 
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for x G [Q,Xm] m+ i], J £i 9 {s,t)ds 

_ 3V (i) V m 

and so 

(*,<)</ \ ei —-2-{ a ,t)\da < — B^(s)ds<C^^. 



dx 3 ' Jj, 9a; 4 ' e 

In a similar way, it can be shown that 

\s t ^p(x,t)\<CB^(x). 
The proof for the and their derivatives is similar. ■ 

Lemma 14. Assume that d r > 0. Let A(x) satisfy 0) and {3J). Then, for 
each i, 1 < i < n., and eac/i (xj,tk) € fl M ' N 

\(L M ' N (W^ - w^fo.fr))! < CfM- 1 + iff (s^) + C^'- 1 ] (69) 



|(L^(W L -w^(^ fc ))| < C[M-'+ £ ,V i? ' (l] ' l) +^ Ij+1 r j'- 1 ]. 

(70) 

Assuming, furthermore, that Xj (f. J, then 

{{L M,N {w J,_ w L )i{xj)tk))l < C[M -! + i ff (g ) + ~ ^i-l) 2 ], (71) 

£r+l 

Analogous results hold for the W R — wf and their derivatives. 

Proof. Using (|4"9")l . (j5Tj) and the bound in Lemma [71 for any (xj,tk) € i7 M,Ar , 

|(L M ' iV (W L -w i ) i ( a;j ,t fc ))| < C[(t fc -t fc _ 1 ) + | ei (^-^)«;f (^,t fc )|]. (72) 

From the decompositions and bounds in Lemma [T31 with (|52"j) and (JSHJ) , it 
follows from (172]) that 
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< C[Eg=l h( D l - &H,<z(Zj>*fc)l + \£i{D 2 x - £2)w l:r+1 {x j ,t k )} 

< C[Ej=i max se/j \siW^(s,t k )\ + (x j+ i - i 3 _i)max s6 ^ \£iW^ +1 (s,t k )\)} 

(73) 



Substituting 1 for each of the min expressions gives ([69]) and (|70|) is obtained 
by substituting the appropriate ratio — in each such expression. 

In the remaining case when Xj ^ J, (|54|) can be used instead of (f53|) . and it 
follows by a similar argument to the above that 

h(% ~ %z)*t{xi,t*)\ < C[£minAl}i^-i) + ( ^ +1 ~ Xi - l)2 ]. 

g— 1 y 

(74) 

Substituting 1 for the min expression, as before, gives (fTTj) . 
The proof for the wf- and their derivatives is similar. ■ 

Lemma 15. Let A(x,t) satisfy 0) and 0). Then, on Q KLN , for each 1 < 
i < n, the following estimates hold 

|(l m, W(w l _ w ^ Mxj>tk)l < C(M -i + B^x^O). (75) 
j4n analogous result holds for the wf~. 

Proof. From (f52|) and Lemma for each i = 1 , . . . , n , it follows that on 

a M < N , 

\(L M ' N (W^ -w^Mx j,t k )\ = \((£ t -Dt)-E(-^-6 2 x ))wf( Xj ,t k )\ 

<C{M- 1 +e i Y d B ' l{Xj - l) ) 
< CiM- 1 + Bfai-i)). 
The proof for the and their derivatives is similar. ■ 

The following theorem provides the error estimate for the singular compo- 
nent. 

Theorem 2. Let A(x,t) satisfy (0) and 0. Letw denote the singular com- 
ponent of the exact solution from (OP and W the singular component of the 
discrete solution from fr47\ ). Then 

||W — w|| < C(M _1 + N~ 2 (lnN) s ). (76) 
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Proof. Since w = w L + w R , it suffices to prove the resuit for w L and w R 
separateiy. Here it is proved for w L by an application of Lemma [TT1 A similar 
proof holds for w R . 

The proof is in two parts: Xj is such that either Xj <£ J or Xj = o> G J. 
First assume that Xj ^ J. Each open subinterval (<7fc,<7fc+i) is treated sepa- 
rately. 

First, consider Xj € (0, cxx). Then, on each mesh M, Xj+i — Xj-i < CN~ 1 ai 
and the result follows from (40} and Lemma fT2l 

Secondly, consider xj 6 (<ti, 02), then a\ < Xj-% and Xj+i —Xj-x < CN~ 1 o-2- 
The 2" possible meshes are divided into subclasses of two types. On the 
meshes fi ' with b\ = the result follows from (14"0"1) . (4"2"} and Lemma 

[T2l On the meshes fi M ' N with 61 = 1 the result follows from (4"0"1) . (f^5|) and 
Lemma 1131 

Thirdly, in the general case Xj £ (cr m ,<7 m +i) for 2 < m < n — 1, it follows 

_i — M,N 

that er m < Xj-i and 2^+1 — < CW a m +i- Then J7 is divided into 

— M,N r-p^M,N , , -p^M,N 

subclasses of three types: J/ = {fi : b\ = ■ ■ ■ = b m = 0), !l r = 

{fi : b r = 1, fe r +i = • • • = b m = for some 1 < r < m — 1} and 

= \T2 M,N : b m = 1}. On 7?Q /:Ar the result follows from (4TJ1) . (35} and 

Lemma [T2l on 7J r ,iV from (40}, (42]), 03} and Lemma [TBI on 72 m ' from 
gD}, (43} and Lemma El 

Finally, for Xj <E (cr n , 1), cr„ < and — aij-i < CN^ 1 . Then J7 
is divided into subclasses of three types: fi^ LN — {fi M ' N : 61 = • • • = b n = 
0}, J7 r ' = {fi ' : b r = 1, 6,+i = • • • = 6„ = for some 1 < r < n — 1} 
and 72 n ' = {72 ' : 6 n = 1}. On 7J ' the result follows from (40}, (42} 
and Lemma [T2] on J? r ' from (|4U} . ((31}, (03} and Lemma [TBI on ' from 



(43} and Lemma [151 

Now assume that Xj G J. Then asj = oy, some r. It suffices to show that 
\(L M > N (W L -w L M Xj ,t k )\ 



\(L M ' N $>)i(xj,tk)\ 



< C, (77) 



for each i = 1, ... , n, because an application of the Comparison Principle 
then yields the required result. 

The bounds on the denominator are given in (59}. To bound the numerator 
note that either d r > or d r — 0. 

Suppose first that d r > 0. Then the cases i > r and i < r are treated 
separately. 

If i > r, then, by (69} in Lemma [T4l 



\(L M ' N (W L - w^fo.t*)! < Cpf" 1 + ^_x) + (78) 
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Since d r > 0, by Lemma[l B^(x^x) = B^{a r - h r ) < CB^{a r ) = CN- 2 , 
and so 

\(L M > N (W L w L Mxj,t k )\ < CIA/- 1 + TV- 2 + '±±^±1], (79 ) 

Using (j^Tj) and the ordering of the £,;, these bounds on the numerator and 

denominator lead to (f77|) . 

If i < r, then, by (JTOJ) in Lemma [Til 

\(L M ' N (W L -^ L )U Xj ,t k )\ < C[M-'+e t y / B " {X ^ l) + J ^ Xj+1 ^ 1 }- 

(80) 

Since d r > 0, by Lemma[8l Xj— i = <r r — h r > x q ^ for 1 < < r — 1 and 



Thus 

\(L M ' N (W L -w L )) i (x j ,t k )\ < CIM- 1 + £ -^N- 2 + £l gj±j_Zg=i ], (81) 

£ r £r+l V £ ' r + 1 

Using (|41[) and the ordering of the Ei, these bounds on the numerator and 
denominator lead to (|77|) . 

Now suppose that d r = 0. Then c? r _i > and D + h r < 0, because otherwise 
Xj (£ J. The cases i > r and i < r are now treated separately. 
If i > r, then, by (I69[) in Lemma [HI with r replaced by r — 1, 

| (£ M,iV (w I, _ w L )) , (x , itfc)| < g[M -l +g L_ i(a .._ l) + ^+1 -^-1 ] ( 82) 



Since d r _i > 0, by Lemma [8] with r replaced by r — 1, 

Br-i( x j-i) = B^_ x {a r - h r ) < CBL_ x (a r -i) = CN~ 2 , 

and so 

\(L M < N (W L - w L Mxj,t k )\ < CiM- 1 + N- 2 + hr + / i r+1 ]. (83) 



Using (|41[) and the ordering of the £j, these bounds on the numerator and 
denominator lead to ([77)1 . 

If z < r, then by (|70[) in Lemma [14] with r replaced by r — 1, 
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\(L M > N (W L w L M Xj ,t k )\ < C[M^+ £l 2 + £i ^+i-^-i ]- 

(84) 

Since d r _i > 0, by Lemma [SJvith r replaced by r — 1, 
for 1 < q < r — 1 and 

£g £V — 1 £V — 1 £y_ 1 

Thus 

| (i M,w (W L_ w L ))t(a ,^ fc) | < C [ M -i + ^_jv-2 + £i^±l^tl]. ( 85 ) 

Using (pLlj) and the ordering of the these bounds on the numerator and 
denominator lead to (l77|) . This completes the proof. ■ 

The following theorem gives the required first order in time and essentially 
second order in space parameter-uniform error estimate. 

Theorem 3. Let A(x,t) satisfy J]|) and (0). Let u denote the exact solution 
of rip and U the discrete solution of J^7[ ). Then 

||U - u|| < CN~ 2 (lnN) 3 . (86) 

Proof. An application of the triangle inequality and the results of Theorems 
Q] and [2] leads immediately to the required result. ■ 
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